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5^ , Abstract 

Motivated by the 'i/Nc expansion, we present a simple model of vrvr scattering as 
a sum of a current- algebra contact term and resonant pole exchanges. The model 
preserves crossing symmetry as well as unitarity up to 1.2 GeV . Key features include 
chiral dynamics, vector meson dominance, a broad low energy scalar (cr) meson and 
a Ramsauer- Townsend mechanism for the understanding of the 980 MeV region. We 
discuss in detail the regularization (corresponding to rescattering effects) necessary to 
make all these nice features work. 
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1 Introduction 



Historically, the analysis of tttt scattering has been considered an important test of our under- 
standing of strong interaction physics (QCD, now) at low energies. It is commonly accepted 
that the key feature is the approximate spontaneous breaking of chiral symmetry. Of course, 
the kinematical requirements of unitarity and crossing symmetry should be respected. The 
chiral perturbation scheme which improves the tree Lagrangian approach by including 
loop corrections and counterterms, can provide a description of the scattering up to the 
energy region slightly above threshold (400 — 500 MeV). 

In order to describe the scattering up to energies beyond this region (say to around 
1 GeV) it is clear that the effects of particles lying in this region must be included and some 
new principle invoked. A plausible hint comes from the large Nc approximation to QCD, 
in which the leading order scattering amplitudes consist of just tree diagrams containing 
resonance exchanges as well as possible contact diagrams 0. The method suggests that an 
infinite number of resonances are required and also a connection with some kind of string 
theory 0. 

Some encouraging features were previously found in an approach which truncated the 
particles appearing in the effective Lagrangian to those with masses up to an energy slightly 
greater than the range of interest. This seems reasonable phenomenologically and is what 
one usually does in setting up an effective Lagrangian. The most famous example is the 
chiral Lagrangian of only pions. In Ref. this Lagrangian provided, as a starting point, 
a contact term which described the threshold region. However the usual observation was 
made that the real part of the / = 0, J = partial wave amplitude quite soon violated the 
unitarity bound \Rq\ < 1/2 rather severely. The inclusion of the contribution coming from 
the p meson exchange was observed to greatly improve, although not completely solve, this 
problem. These results are shown explicitly in Fig. 1 and provide some encouragement for 
the possible success of a truncation scheme. 

In Ref. [^, it was observed that the inclusion of resonances up till and including the 
p-wave region enabled one to construct an amplitude which satisfied the unitarity bounds 
up to about L3 GeV. It was assumed that, above this point, new resonances would come 
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Figure 1: Predicted curves for Rq. The solid line which shows the current algebra 
+ p result for is much closer to the unitarity bound of 0.5 than the dashed 
line which shows the current algebra result alone. 
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in to preserve unitarity. This hypothesis was called local cancellation. The model produced 
a reasonable looking I = J = phase shift up to about 800 MeV. In this paper we will 
attempt to describe and carefully compare with experiment the interesting physics lying 
between 800 and 1200 MeV in this truncated 1/Nc inspired framework. Specifically we will 
be concerned with the proper inclusion of the /o(980) scalar resonance as well as the opening 
of the KK channel. We find that a simple reasonable description of the /o(980) region is 
obtained when the interplay of this resonance with its background is taken into account. In 
this approach the background amplitude is predicted by the model itself. In the region just 
above the KK threshold we notice the feature analogous to the elastic case that the severe 
unitarity violation of the inelastic tttt — > KK amplitude is damped by the inclusion of vector 
meson and scalar meson exchange diagrams. 

Of course, it would be wonderful if one could simply add the various contributions to the 
tree level amplitude and find a good match to experiment. This is not possible for a variety 
of reasons, which are discussed in Section 2. The needed regularizations are introduced there. 
Section 3 gives a brief overview of the model and reviews the important role of a broad scalar 
meson in the low energy (< 800 MeV) region. Section 4 contains a discussion of various 
aspects of the 1 GeV region. The characteristic feature - a type of Rams auer- Towns end 
effect resulting from the interplay of the /o(980) resonance with the predicted background 
- is outlined in section 4.1 and treated in more detail in 4.2. In section 4.3 it is shown 
that the introduction of the next group of resonances, located in the 1300 MeV region, does 
not make major changes in the tttt scattering below 1200 MeV (the changes are essentially 
absorbed in small changes of the parameters of the broad low energy scalar). In section 
4.4 it is demonstrated that the phenomenological introduction of inelastic effects associated 
with the opening of the KK channel does not make a significant change in our picture of 
TTTT TTTT bclow 1200 MeV. Section 4.5 contains a presentation of the I = J = phase 
shift obtained by combining our predicted real part with unitarity. In section 5 we discuss 
the inelastic tttt — > KK channel and show that here also the resonance exchanges damp 
the unitarity bound violation due to the contact term. Section 6 contains the summary 
and further discussion. Finally, Appendices A, B and C give details on, respectively, the 
scattering kinematics, the chiral Lagrangian and the unregularized amplitudes. 
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2 Difficulties of the Approach 

In the large picture the leading amplitude (of order 1/Nc) is a sum of polynomial contact 
terms and tree type resonance exchanges. Furthermore the resonances should be of the simple 
qq type; glueball and multi-quark meson resonances are suppressed. In our phenomenological 
model there is no way of knowing a priori whether a given experimental state is actually of 
qq type. For definiteness we will keep all relevant resonances even though the status of a 
low lying scalar resonance like the /o(980) has been considered especially controversial |Q. If 
such resonances turn out in the future to be not of qq type, their tree contributions would be 
of higher order than l/N^. In this event the amplitude would still of course satisfy crossing 
symmetry. 

The most problematic feature involved in comparing the leading 1/Nc amplitude with 
experiment is that it does not satisfy unitarity. In fact, resonance poles like 

Jih <^-^) 

will yield a purely real amplitude, except at the singularity, where they will diverge and 
drastically violate the unitarity bound. Thus in order to compare the 1/Nc amplitude with 
experiment we must regularize the denominators in some way. The usual method, as em- 
ployed in Ref . , is to regularize the propagator so that the resulting partial wave amphtude 
has the locally unitary form 

(2.2) 



M2 - s - iMT ■ 

This is only valid for a narrow resonance in a region where the background is negligible. Note 
that the —iMT is strictly speaking a higher order in 1/Nc effect. 

For a very broad resonance there is no guarantee that such a form is correct. Actually, in 
Ref. 1^ it was found necessary to include a rather broad low lying scalar resonance (denoted 
cr(550)) to avoid violating the unitarity bound. A suitable form turned out to be of the type 

MG 



M'^-s- iMG' 



(2.3) 



where G is not equal to the parameter G' which was introduced to regularize the propagator. 
Here G is the quantity related to the squared coupling constant. 
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Even if the resonance is narrow, the effect of the background may be rather important. 
This seems to be true for the case of the /o(980). Demanding local unitarity in this case 
yields a partial wave amplitude of the well known form |^: 

/ ^^^^ + e'^in5, (2.4) 

where 5 is a background phase (assumed to be slowly varying). We will adopt a point of 
view in which this form is regarded as a kind of regularization of our model. Of course, non 
zero 6 represents a rescattering effect which is of higher order in l/Nc- The quantity e^**^, 
taking 6 = constant, can be incorporated into the squared coupling constant connecting the 
resonance to two pions. In this way, crossing symmetry can be preserved. From its origin, 
it is clear that the complex residue does not signify the existence of a ghost particle. The 
non-pole background term in eq. ( |2.4| ) and hence 6 is to be predicted by the other pieces in 
the effective Lagrangian. 

Another point which must be addressed in comparing the leading l/N^ amplitude with 
experiment is that it is purely real away from the singularities. The regularizations mentioned 
above do introduce some imaginary pieces but these are clearly more model dependent. Thus 
it seems reasonable to compare the real part of our predicted amplitude with the real part 
of the experimental amplitude. Note that the difficulties mentioned above arise only for the 
direct channel poles; the crossed channel poles and contact terms will give purely real finite 
contributions. 

It should be noted that if we predict the real part of the amplitude, the imaginary part 
can always be recovered by assuming elastic unitarity (which is likely to be a reasonable 
approximation up to about 1 GeV) . Specializing eq. ( |A.6[ ) in Appendix A to the tttt channel 



we have for the imaginary piece // of the J, I partial wave amplitude 

1 



1 ± Vr//' - 4R{' 



where 7]^ is the elasticity parameter. Obviously this formula is only meaningful if the real 
part obeys the bound 

\Ri\ < y • (2.6) 
The main difficulty one has to overcome in obtaining a unitary amplitude by the present 
method is the satisfaction of this bound. Therefore, one sees that making regularizations like 







M{MeV) 


TtotiMeV) 


Br{2'n)% 


cr(550) 


0+(0++) 


559 


370 


— 


p(770) 


1+(1-) 


769.9 


151.2 


100 


/o(980) 


0+(0++) 


980 


40-400 


78.1 


/2(1270) 


0+(2++) 


1275 


185 


84.9 


/o(1300) 


0+(0++) 


1000-1500 


150-400 


93.6 


p(1450) 


1+(1-) 


1465 


310 


seen 



Table 1: Resonances included in the tttt tttt channel as listed in the PDG. 
Note that the a is not present in the PDG and is not being described exactly 
as a Breit- Wigner shape; we listed the fitted parameters shown in column 1 of 
Table 2 where G' is the analog of the Breit- Wigner width. 



eqs. ( |2.2|) and ( |2.4| ) which provide unitarity in the immediate region of a narrow resonance 
is not at all tantamount to unitarizing the model by hand. One might glance again at Fig. 1 
for emphasis of this point. 

To summarize this discussion, we will proceed by comparing the real part of a suitably 
regularized tree amplitude computed from a chiral Lagrangian of pseudoscalar mesons and 
resonances with the real part of the experimental amplitude deduced from the standard 
phase shift analysis. 



3 Overview and Low Energy Region 

The amplitude will be constructed from the non-linear chiral Lagrangian briefly summarized 
in Appendix B. To start with, we shall neglect the existence of the K mesons. Then the 
form of the unregularized amplitude is identical to the one presented in Ref. [Q. The neutral 
resonances which can contribute have the quantum numbers J^*" = 0"'"+, 1 , and 2"*"+. We 
show in Table 1 the specific ones which are included, together with their masses and widths, 
when available from the Particle Data Group (PDG) listings. 

Essentially there are only three arbitrary parameters in the whole model, these correspond 
to the three unknowns in the description of a broad scalar resonance given by eq. ( p.3|) . We 
will include only the minimal two derivative chiral contact interaction contained in eq. ( |B.7|) 
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of Appendix B. Clearly, higher derivative contact interaction may also be included (see, for 
example, sec. III.E of Ref. IQ). 

As shown in Fig. 1, although the introduction of the p dramatically improves unitarity 
up to about 2 GeV , violates unitarity to a lesser extent starting around 500 MeV . (As 
noted in Ref. |Q, the J = J = channel is the only troublesome one.) To completely restore 
unitarity in the present framework it is necessary to include a low mass broad scalar state 
which has historically been denoted as the a. It seems helpful to recall the contribution of 
such a particle to the real part of the amplitude component A{s,t,u) defined in eq. ( |A.8| ): 

4 / . N 327r G , ^ {Ml - s) + iM^G' 



where 



2 \ o / 2 \ 2 

ml \ ^ / m; ^ 



and G is related to the coupling constant 70 defined in eq. ( p.ll| ) by 

G = 7o-^- (3.3) 

Note that the factor (s — 2m^)^ is due to the derivative-type coupling required for chiral 
symmetry in eq. ( p.ll|) . The total amplitude will be crossing symmetric since A(s,t,M) 
and y4(M, t,s) in eq. ( |A.8| ) are obtained by performing the indicated permutations. G' is a 
parameter which we introduce to regularize the propagator. It can be called a width, but it 
turns out to be rather large so that, after the p and vr contributions are taken into account, 
the partial wave amplitude 7?° does not clearly display the characteristic resonant behavior. 
In the most general situation one might imagine that G could become complex as in eq. ( |2.4|) 
due to higher order in IjNf, corrections. It should be noted, however, that eq. ( p.4| ) expresses 
nothing more than the assumption of unitarity for a narrow resonance and hence should not 
really be applied to the present broad case. A reasonable fit was found in Ref. |Q for G 
purely real, but not equal to G' . By the use of eq. (^.51), unitarity is in fact locally satisfied. 

A best overall fit is obtained with the parameter choices; = 559 MeV ^ G/G' = 0.29 
and G' = 370 MeV . These have been slightly fine-tuned from the values in Ref. [Q in order 
to obtain a better fit in the 1 GeV region. The result for the real part -Rq due to the inclusion 
of the a contribution along with the vr and p contributions is shown in Fig. 2. It is seen that 
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Figure 2: The solid line is the current algebra + p + a result for Rq. The 
experimental points, in this and suceeding figures, are extracted from the phase 
shifts using eq. ( [A. 61 ) and actually correspond to Rq/tIq. (□) are extracted from 
the data of Ref. [§] while (A) are extracted from the data of Ref. The 
predicted Rq is small around the 1 GeV region. 
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the unitarity bound is satisfied and there is a reasonable agreement with the experimental 
points ^ up to about 800 MeV. Beyond this point the effects of other resonances (mainly 
the /o(980)) are required. From eqs. O and ( [ATTD we see that the contribution 

of a to Rq turns negative when s > M^. This is the mechanism which leads to satisfaction 
of the unitarity bound (c.f. Fig. 1). For s < one gets a positive contribution to Rq. 
This is helpful to push the predicted curve upwards and closer to the experimental results in 
this region, as shown in Fig. 3. The four-derivative contribution in the chiral perturbation 
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Figure 3: A blowup of the low energy region. The solid line is the current 
algebra + p contribution to -Rq. The dashed line includes the a and has the effect 
of turning the curve down to avoid unitarity violation while boosting it at lower 
energies. 



theory approach performs the same function; however it does not change sign and hence 
does not satisfy the unitarity bound above the 450 MeV region . 
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4 The 1 GeV Region 



4.1 The main point 

Reference to Fig. 2 shows that the experimental data for he considerably lower than the 
vr + p + (T contribution between 0.9 and 1.0 GeV and then quickly reverse sign above this 
point. We will now see that this distinctive shape is almost completely explained by the 
inclusion of the relatively narrow scalar resonance /o(980) in a suitable manner. One can 
understand what is going on very simply by starting from the real part of eq. 



^ (M^ - s) cos(26) - MT sm(26) 1 , , , 

This expresses nothing more than the restriction of local unitarity in the case of a narrow 
resonance in the presence of a background. We have seen that the difficulty of comparing the 
tree level 1/Nc amplitude to experiment is enhanced in the neighborhood of a direct channel 
pole. Hence it is probably most reliable to identify the background term - sin(25) with our 
prediction for Rq. In the region of interest, Fig. 2 shows that Rq is very small so that one 
expects, 6 to be roughly 90° (assuming a monotonically increasing phase shift). Hence the 
first, pole term is approximately 

(M2 - sy + M2r2 ' ^ ■ ^ 

which contains a crucial reversal of sign compared to the real part of eq. ( |2.2| ). Thus, just 
below the resonance there is a sudden negative contribution which jumps to a positive one 
above the resonance. This is clearly exactly what is needed to bring experiment and theory 
into agreement up till about 1.2 GeV, as is shown in Fig. 4. The actual amplitude used for 
this calculation properly contains the effects of the pions' derivative coupling to the /o(980) 



as in eq. (|3?T 



It is interesting to contrast this picture with Fig. 10 in Ref. There the interaction 
with the background was not taken into account and there was no reversal of sign. Thus, 
although the unitarity bound was obeyed, the experimental phase shifts could only be prop- 
erly predicted up to about 0.8 GeV. If the /o(980) contribution in that Fig. 10 is flipped in 
sign it is seen to agree with the present Fig. 4. 
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Figure 4: (a): The solid line is the current algebra + p + a + /o(980) result 
for Rq obtained by assuming column 1 in Table 2 for the a and /o(980) param- 
eters (i?r(/o(980) — > 27t) = 100%). (b): The solid line is the current algebra 
+ p + a + /o(980) result for Rq obtained by assuming column 2 in Table 2 
(Sr(/o(980) ^ 27i) = 78.1%) . 



The above mechanism, which leads to a sharp dip in the I = J = partial wave con- 
tribution to the TTTT-scattering cross section, can be identified with the very old Ramsauer- 



Townsend effect [|TT| which concerned the scattering of 0.7 eV electrons on rare gas atoms. 
The dip occurs because the background phase of 7r/2 causes the phase shift to go through vr 
(rather than tt/2) at the resonance position. (Of, course, the cross section is proportional to 
J2i j(2t/+ 1) sin^(5/).) This simple mechanism seems to be all that is required to understand 
the main feature of tttt scattering in the 1 GeV region. 

4.2 Detailed analysis 

Here we will compare with experimental data, the real part of the I = J = partial wave 
amplitude which results from our crossing symmetric model. First we will consider the sum 
of the contributions of the current algebra, p-meson, a and /o(980) pieces. Then we will add 
pieces corresponding to the next group of resonances; namely, the /2(1270), the p{1450) and 
the /o(1300). In this section we will continue to neglect the KK channel. 
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The current algebra plus p contribution to the quantity A{s, t, u) defined in eq. ( |A.8 



F2 



y pTTTV 



U 



[m? — t) — inipT p9{t — Ami] 



+ 



{m?p — u) — impV p6{u — Ami) 



(4.3) 



Note that for the / = J = channel this will yield a purely real contribution to the partial 
wave amplitude. The contribution of the low lying a meson was given in eq. ( |3.1| ). For the 
important /o(980) piece we have 



ReA 



/o(980) 



[s, t, u) = Re 



2i5/ 



2m. 



2\2 



mj^- s- imf^TtotifoWis - Ami) 



(4.4) 



where 5 is a background phase parameter and the real coupling constant •yfgnn is related to 
the /o(980) nn width by 



r(/o(980) ^ tttt) 



3 If 
32tt mjg \^ 



Ami 



m 



(4.5) 



/o 



We will not consider 5 to be a new parameter but shall predict it as 

^sm{25) ^ R',{s = ml) , 



(4.6) 



where Rq is computed as the sum of the current algebra, p, and sigma pieces. Since the 
KK channel is being neglected, one might want to set the regularization parameter Ttotifo) 
in the denominator to r(/o(980) tttt). We shall try both this possibility as well as the 



experimental one 



r(/o(980) 



TTTT 



78.1%. 



^totifo) 

A best fit of our parameters to the experimental data results in the curves shown in 
Fig. 4 for both choices of branching ratio. Only the three parameters G/G', G' and Mg- 
are essentially free. The others are restricted by experiment. Unfortunately the total width 
Ttotifo) has a large uncertainty; it is claimed by the PDG to lie in the 40 — 400 MeV range. 
Hence this is effectively a new parameter. In addition we have considered the precise value 
of m/o to be a parameter for fitting purposes. The parameter values for each fit are given in 



§We introduced the step function 9{s — Aral) ™ the propagator and have checked that its inclusion does 
not make much difference in the results. 
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With Next Group 




No p(1450) 


Si?(/o(980) ^ 27r)% 


100 


78.1 


78.1 


78.1 


100 


78.1 


78.1 


78.1 


100 


%° 


1 


1 


0.8 


0.6 


1 


1 


0.8 


0.6 


1 


Mj„(98o) {MeV) 


987 


989 


990 


993 


991 


992 


993 


998 


992 


Ttot {MeV) 


64.6 


77.1 


75.9 


76.8 


66.7 


77.2 


78.0 


84.0 


64.6 


Mo- (Me\/) 


559 


557 


557 


556 


537 


537 


535 


o o 

533 


525 


G" {MeV) 


370 


371 


380 


395 


422 


412 


426 


451 


467 


G/G' 


0.290 


0.294 


0.294 


0.294 


0.270 


0.277 


0.275 


0.270 


0.263 


5 (deg.) 


85.2 


86.4 


87.6 


89.6 


89.2 


89.7 


91.3 


94.4 


90.4 




2.0 


2.8 


2.7 


3.1 


2.4 


3.2 


3.2 


3.4 


2.5 



Table 2: Fitted parameters for different cases of interest. 

Table 2 together with the values. It is clear that the fits are good and that the parameters 
are stable against variation of the branching ratio. The predicted background phase is seen 
to be close to 90° in both cases. Note that the fitted width of the /o(980) is near the low 
end of the experimental range. The low lying sigma has a mass of around 560 MeV and a 
width of about 370 MeV . As explained in section 3, we are not using exactly a conventional 
Breit- Wigner type form for this very broad resonance. The numbers characterizing it do 
however seem reasonably consistent with other determinations [12], |T^. 



4.3 Effect of the next group of resonances 

Going up in energy we encounter J^^ = 2"'"+,0^"'" and 1 resonances in the 1300 MeV 
region. The properties of the 2^~^ state /2(1270) are very well established. For the others 
there is more uncertainty but the PDG lists the /o(1300) and p(1450) as established states. 
However the mass of the /o(1300) can apparently lie anywhere in the 1000 — 1500 MeV 
range. In Ref. @] it was noted that the contributions of these next group particles tended to 
cancel among themselves. Thus we do not expect their inclusion to significantly change the 
previous results in the range of interest up to about 1.2 GeV. 

In Fig. 5 we display the contribution of the next group particles by themselves to Rq. (The 
amplitudes are summarized in Appendix C). The dashed curve is essentially a reproduction 
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Figure 5: Contribution from the next group of resonances; the sohd curve is 
obtained with the reverse sign of the /o(1300) piece. 

of Fig. 6 of Ref. [^]. The somewhat positive net contribution of these resonances to -Rg is 
compensated by readjustment of the parameters describing the low lying sigma. It may be 
interesting to include the effect of the background phase for the /o(1300) as we have just 
seen that it was very important for the proper understanding of the /o(980). To test this 
possibility we reversed the sign of the /o(1300) contribution and show the result as the solid 
curve in Fig. 5. This sign reversal is reasonable since our model suggests a background phase 
of about 270° in the vicinity of the /o(1300). It can be seen that there is now a significantly 
greater cancellation of the next group particles among themselves up to about 1.2 GeV. The 
resulting total fits are shown in Fig. 6 for both 100% and 78.1% assumed /o(980) vrvr 
branching ratios and the parameters associated with the fits are shown in Table 2. It is clear 
that the fitted parameters and results up to about 1.2 GeV are very similar to the cases when 
the next group was absent. Above this region, there is now, however, a positive bump in Rq 
at around 1.3 GeV. This could be pushed further up by choosing a higher mass (within the 
allowable experimental range) for the /o(1300). Resonances in the 1500 MeV region, which 
have not been taken into account here, would presumably also have an important effect in 
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Ys (GeV) NTs" (GeV) 
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Figure 6: Prediction for i?Q with the next group of resonances, (a) assumes 
(column 5 in Table 2) {B R{f oi980) 27r) = 100%) while (b) assumes (column 
6) (5i?(/o(980) ^ 27r) = 78.1%). 

the region above 1.2 GeV. Clearly, there is not much sense, at the present stage, in trying 
to produce a fit above 1.2 GeV. 

The analysis above assumed that the p(1450) decays predominantly into two pions since 
the PDG listing does not give any specific numbers. On the other hand the K*[1A10), which 
presumably is in the same SU (3) multiplet as the p(1450), has only a 7% branching ratio into 
Kn. Thus it is possible that p(1450) actually has a small couphng to tttt. To test this out we 
redid the calculation with the complete neglect of the p(1450) contribution. The resulting fit 
is shown in the last column of Table 2 and it is seen to leave the other parameters essentially 
unchanged. 

It thus seems that the results are consistent with the hypothesis of local cancellation, 
wherein the physics up to a certain energy E is described by including only those resonances 
up to slightly more than E and it is furthermore hypothesized that the individual particles 
cancel in such a way that unitarity is maintained. 
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4.4 Effects of inelasticity 

Up to now we have completely neglected the effects of coupled inelastic channels. Of course 
the 47r channel opens at 540 MeV, the Gtt channel opens at 810 MeV and, probably most 
significantly, the KK channel opens at 990 MeV. We have seen that a nice undestanding of 
the TT7C elastic channel up to about 1.2 GeV can be gotten with complete disregard of inelastic 
effects. Nevertheless it is interesting to see how our results would change if experimental data 
on the elasticity parameter tjq are folded into the analysis. Figure [7| illustrates the results for 



0.8 



0.6 




0.4 - 



0.2 



0.0 



Figure 7: An experimental determination of ?7q = a/1 — 4|T{2^oP Ell- 



77q(s) obtained from an experimental analysis [|1^ of vrvr — » KK scattering. For simplicity, we 
approximated the data by a constant value r/g = 0.8 above the KK threshold. Figure 8(a) 
shows the effect of this choice on -Ro(s) computed without the inclusion of the next group 
of resonances, while Fig. 8(b) shows the effect when the next group is included. Comparing 
with Fig. 4(b) and 6(b), we see that setting rj^ = 0.8 has not made any substantial change. 
The parameters of the fit are shown in Table 2 as are the parameters for an alternative fit 
with ?7q = 0.6. The latter choice leads to a worse fit for Rq. 

We conclude that inelastic effects are not very important for understanding the main 
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(a) (b) 

Figure 8: Predictions with phenomenological treatment of inelasticity (r^Q = 0.8) 
above KK threshold, (a): without next group, (b): with next group. 



features of tttt scattering up to about 1.2 GeV . However, we will discuss the calculation of 
?7q(s) from our model in section 5. 

4.5 Phase shift 

Strictly speaking our initial assumption only entitles us to compare, as we have already done, 
the real part of the predicted amplitude with the real part of the amplitude deduced from 
experiment. Since the predicted -Ro(s) up to 1.2 GeV satisfies the unitarity bound (within 
the fitting error) we can calculate the imaginary part Iq{s), and hence the phase shift 5q{s) 
on the assumption that full unitarity holds. This is implemented by substituting R^{s) into 
eq. (p75|) and resolving the discrete sign ambiguities by demanding that 5q{s) be continuous 
and monotonically increasing (to agree with experiment). It is also necessary to know ?7o(s) 
for this purpose; we will be content with the approximations above which seem sufficient for 
understanding the main features of vrvr scattering up to 1.2 GeV . 

In this procedure there is a practical subtlety already discussed at the end of section 
IV of Ref. [0. In order for 6q{s) to increase monotonically it is necessary that the sign in 
front of the square root in eq. ( p3| ) change. This can lead to a discontinuity unless 2|_Rq(s)| 
precisely reaches i]q{s). However the phase shift is rather sensitive to small deviations from 
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this exact matching. Since the fitting procedure does not enforce that |-Ro(s)| go precisely to 
?7q(s)/2 pa 0.5, this results in some small discontinuities. (These could be avoided by trying 
to fit the phase shift directly.) 

Figure 9 shows the phase shift Sq{s) estimated in this manner for parameters in the first 




fV (GeV) 

Figure 9: Estimated phase shift using the predicted real part and unitarity rela- 
tion. 

column of Table 2. As expected, the agreement is reasonable. A very similar estimate is 
obtained when (column 3 of Table 2) tJq is taken to be 0.8 while considering the tttt branching 
ratio of /o(980) to be its experimental value of 78.1%. It appears that these two parameter 
changes are compensating each other so that one may again conclude that the turning on 
of the KK channel really does not have a major effect. When the next group of resonances 
is included (column 7 of Table 2) the estimated Sq{s) is very similar up to about 1.2 GeV. 
Beyond this point it is actually somewhat worse, as we would expect by comparing Fig. 8(b) 
with Fig. 8(a). 
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5 7T7T KK Channel 

We have seen that tttt scattering can be understood up to about 1.2 GeV with 

the neglect of this inelastic channel. In particular, a phenomenological description of the 
inelasticity did not change the overall picture. However we would like to begin to explore 
the predictions of the present model for this channel also. The whole coupled channel 
problem is a very complicated one so we will be satisfied here to check that the procedure 
followed for the vrvr elastic channel can lead to an inelastic amplitude which also satisfies the 
unitarity bounds. Specifically we will confine our attention to the real part of the I = J = 
TTTT KK amplitude, -R?2 o defined in eq. (|A.11|) . 

In exact analogy to the tttt tttt case we first consider the contribution of the contact 
plus the i^*(892) plus the o"(550) terms. It is necessary to know the coupling strength of the 
a to KK, defined by the effective Lagrangian piece 

-l^ad^Kd^K. (5.1) 

If the a is ideally mixed and there is no OZI rule violating piece we would have 'j^^^j^ = 7o 
as defined in eq. ( |B.11| ). For definiteness, we shall adopt this standard mixing assumption. 
The appropriate amplitudes are listed in Appendix C. Figure 10 shows the plots of -R?2;0 
the current algebra part alone, the current algebra plus K* and the current algebra plus K* 
plus cr parts. Notice that unitarity requires 




The current algebra result already clearly violates this bound at 1.05 GeV. As before, this is 
improved by the K* vector meson exchange contribution and further improved by the very 
important tail of the a contribution. The sum of all three shows a structure similar to the 
corresponding Fig. 2 in the tttt —>■ tttc case. The unitarity bound is not violated until about 
1.55 GeV. 

Next, let us consider the contribution of the /o(980) which, since the resonance straddles 
the threshold, is expected to be important. We need to know the effective coupling constant 
of the /o to TTTT and to KK. As we saw in eq. ( ^4.4| ), and the subsequent discussion, the 
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Figure 10: Contributions to tttt — > KK (-R?2;o)- solid line shows the current 
algebra result; the dashed line represents the inclusion of K*{S92); the dotted 
line includes the (t(550) too. 
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effective vrvr coupling should be taken as 7/o7r7re*^. Experimentally, only the branching ratios 

for /o(980) nn and /o(980) KK are accurately known. We will adopt for definiteness 

the value of 7/o7r7r corresponding to the fit in the third column of Table 2 (rtot(/o(980)) = 

76 MeV). It is more difficult to estimate the /o(980) — > KK effective coupling constant since 

the central value of the resonance may actually lie below the threshold. By taking account i 

of the finite width of the /o(980) we get the rough estimate |7jp^;^| = 10 GeV~^ ~ 4|7/(,^7r| 

for the choice in the third column, M/q(98o) = 990 MeV. Of course, this estimate is very 

sensitive to the exact value used for Mjj,(98o). It seems reasonable to take Ij^kk t)e purely 

real. The results of including the /o(980) contribution, for both sign choices of 'Jj^xk^ 

shown in Fig. 11. The unitarity bounds are satisfied for the positive sign of If^^xT? ^lut 

slightly violated for the negative sign choice. 

Finally, let us consider the contributions to vrvr — ^ KK from the members of the multiplets 

containing the next group of particles. There will be a crossed channel contribution from 

the strange excited vector meson i^'*(1410). However it will be very small since i^'*(1410) 

predominately couples to K*tt and has only a 7% branching ratio to i^vr. In addition there 

will be a crossed channel scalar irQ(1430) diagram as well as a direct channel scalar /o(1300) 

diagram contributing to tttt KK. The /o(1300) piece is small because /o(1300) has a very 

small branching ratio to KK. Furthermore the ii'o(1430) piece turns out also to be small; 

we have seen that the crossed channel scalar gave a negligible contribution to tttt tttt. 

The dominant next group diagrams involve the tensor mesons. Near threshold, the crossed 

channel 7^^2(1430) diagram is the essential one since the direct channel /2(1270) contribution 

for the J = partial wave is suppressed by a spin-2 projection operator. Above 1270 MeV 

^With rtot(/o(980)) = 76 MeV we would have r(/o(980) -> KK)) = 16.6 MeV. Then 7^^^^ is estimated 
from the formula: 

16.6 MeV = |7/„/^kI' / p{M)\A{h{M) -> KK)\''^{M) dM , 

where A(/o(Af) KK) is the reduced amplitude for an /o of mass M to decay to KK, $(M) is the phase 
space factor and p{M) is the weighting function given by 

p{M) = J--^exp\~2 

V TT 1 tot L 

Here, Mq is the central mass value of the /o(980). 



(M - Mo) 



2 



tot 
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the /2(1270) contribution becomes increasingly important although it has the opposite sign to 
the crossed channel tensor piece. Figure 12 shows the net prediction for -R?2 o obtained with 
the inclusion of the main next group contributions from the i^'2(1430) and /2(1270). Both 
assumed signs for If^^li shown and other parameters correspond to column 3 of Table 
2. Clearly there is an appreciable effect. Figure 13 shows the magnitude of |-R?2:ol together 



with one experimental determination [jT4[ of \T^2-fi\ = y (-^i2;o)^ + (-^i2;o)^- The positive sign 
of IfgK'K is favored but, considering the uncertainty in {'Jj^^xkI among other things, we shall 
not insist on this. It seems to us that the main conclusion is that the unitarity bound can 
be satisfied in the energy range of interest. 



6 Summary and Discussion 

We have obtained a simple approximate analytic form for the real part of the tttt scattering 
amplitude in the energy range from threshold to about 1.2 GeV . It satisfies both crossing 
symmetry and (more non-trivially) unitarity in this range. Inspired by the leading 1/Nc 
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Figure 12: Effects on tttt KK due to the next group of resonances for tlie two 
different sign clioices in Fig. 11. 
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approximation, we have written the amphtude as the sum of a contact term and poles. 
Of course the leading 1/Nc amplitude can not be directly compared with experiment since 
it is purely real (away from the direct channel poles) and diverges at the pole positions. 
Furthermore, an infinite number of poles, and higher derivative interactions are in principle 
needed. To overcome these problems we have employed the following procedure. 

a. We specialized to predicting the real part of the amplitude. 

b. We postulated that including only resonances from threshold to slightly more than 
the maximum energy of interest is sufficient. We have seen that this local cancellation 
appears stable under the addition of resonances in the 1300 MeV range. Beyond this 
range we would expect still higher resonances to add in such a way so as to enforce 
unitarity at still higher energies. 

c. In the effective interaction Lagrangian we included only terms with the minimal number 
of derivatives consistent with the assumed chiral symmetry. 

d. The most subtle aspect concerns the method for regularizing the divergences at the 
direct channel resonance poles. In the simplest case of a single resonance dominating 
a particular channel (e.g. the p meson) it is sufficient to add the standard width term 
to the denominator (e.g. the real part of eq. ( p.2|) ). For an extremely broad resonance 
(like a needed low energy scalar isosinglet) the concept of width is not so clear and 
we employed the slight modification of the Breit-Wigner amplitude given in eq. ( |2.3| ). 
Finally, for a relatively narrow resonance in the presence of a non-negligible background 
we employed the regularization given in eq. ( p.4| ) which includes the background phase. 
Self-consistency is assured by requiring that the background phase should be predicted 
by the model itself. 

All the regularizations introduced above are formally of higher than leading order in the 
1/Nc expansion (i.e. of order l/N^ and higher) and correspond physically to rescattering 
effects. In the case of non-negligible background phase, there is an interesting difference from 
the usual tree- level treatment of pole diagrams. The effective squared coupling constant, 
g'piT^T^ of such a resonance to two pions, is then not necessarily real positive. Since this 
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regularization is interpreted as a rescattering effect it does not mean tliat gliost fields are 
present in tlie tlieory. Tliis formulation maintains crossing symmetry wliicli is typically lost 
wlien a unitarization metliod is employed. 

In this analysis, the most non-trivial point is the satisfaction of the unitarity bound for 
the predicted real part of the partial wave amlitude, 

\Ri\<Y, (6.1) 

where 7]^ < 1 is the elasticity parameter. The well known difficulty concerns Rq. If //f (s) is 
known or calculated the imaginary part //(s) can be obtained, up to discrete ambiguities, 
byeq. (|2j). 

The picture of tttt scattering in the threshold to slightly more than 1 GeV range which 
emerges from this model has four parts. Very near threshold the current algebra contact term 
approximates Rq{s) very well. The imaginary part Iq{s), which is formally of order 
can be obtained from unitarity directly using eq. ( p. 51 ) or, equivalently, by chiral perturbtion 
theory. At somewhat higher energies the most prominent feature is the p meson pole in the 
I = J = 1 channel. The crossed channel p exchange is also extremely important in taming 
the elastic unitarity violation associated with the current algebra contact term (Fig. 1). 
Even with the p present. Fig. 1 shows that unitarity is still violated, though much less 
drastically. This problem is overcome by introducing a low mass ~ 550 MeV, extremely 
broad sigma meson. It also has another desirable feature: Rq{s) is boosted (see Fig. 3) closer 
to experiment in the 400 — 500 MeV range. The three parameters characterizing this particle 
are essentially the only unknowns in the model and were determined by making a best fit. 
In the 1 GeV region it seems clear that the /o(980) resonance, interacting with the predicted 
background in the manner of the Ramsauer-Townsend eSect, dominates the structure of the 
I = J = phase shift. The inelasticity associated with the opening of the KK threshold has 
a relatively small effect. However we also presented a preliminary calculation which shows 
that the present approach satisfies the unitarity bounds in the inelastic vrvr KK channel. 

Other recent works [l^, |15|, |16| which approach the problem in different ways, 
also contain a low mass broad sigma. The question of whether the lighter scalar mesons 
are of qq type or meson-meson type has also been discussed P, |12|, |13|. In our model it 
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is difficult to decide this issue. Of course, it is not a clean question from a field theoretic 
standpoint. This question is important for understanding whether the contributions of such 
resonances are formally leading in the 1/Nc expansion. We are postponing the answer as 
well as the answer to how to derive the rescattering effects that were used to regularize the 
amplitude near the direct channel poles as higher order in 1/Nc corrections. Presumably, the 
rescattering effects could some day be calculated as loop corrections with a (very complicated) 
effective Wilsonian action. This would be a generalization of the chiral perturbation scheme 
of pions. Another aspect of the 1/Nc picture concerns the infinite number of resonances 
which are expected to contribute already at leading order. One may hope that the idea of 
local cancellation will help in the development of a simple picture at high energies which 
might get patched together with the present one. Is the simple high energy theory a kind of 
string model ? 

From a practical standpoint (without worrying about all the theoretical issues involved in 
making a comparison with the l/Nc expansion) we have demonstrated that it is possible to 
understand tttt scattering up to the 1 GeV region by shoehorning together poles and contact 
term contributions employing a suitable regularization procedure. It seems likely that any 
crossing symmetric approximation will have a similar form. This is in the spirit of mean field 
theories. 
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Appendix A 

Scattering kinematics 

The general partial wave scattering matrix for the multi channel case can be written as: 



Sab = 5ab + 2lTab . (A.l) 

For simplicity, the diagonal isospin and angular momentum labels have not been indicated. 

By requiring the unitarity condition S"^ S = 1 one deduces for the two channel case the 
following relations: 



/m(Tn) 


= \Tn 




\T2l\ 




= \T22 




\Tl2\ 



(A.2) 

Im^Tu) = T^i Ti2 + Tj'sTss , 

where T12 = T21. In the present case we will identify 1 as the tttt channel and 2 as the KK 
channel. In order to get the relations between the relative phase shifts and the amplitude 
we need to consider the following parameterization of the scattering amplitude: 

S = ' ^ ^ ' , (A.3) 



where S^^k = + and < 77 < 1 is the elasticity parameter. By comparing eq. ( |A.3| ) and 
eq. ( [A.l| ) one can easily deduce: 



r]^ = 1 - 4\T,2\^ . (A.4) 

Analogously, for Taa we have: 

TLM = SUlS^^ ^ , (A.5) 

where I and / label the angular momentum and isospin, respectively. Extracting the real 
and imaginary parts via 

1 - r?/ cos(2(5^) 

IL, = ' ^ ^ "''^ (A.6) 
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leads to the very important bounds 

K;/l<^, 0</i^,<l. (A.7) 

The unitarity also requires iT/g.J < 1/2 . 

Now we relate these partial wave amplitudes to the invariant amplitudes. The invariant 
amplitude for Tii^pi) + T^j{p2) T^kipz) + T^iiPi) is decomposed as: 

5ij5kiA{s, t, u) + 5ik5jiA{t, s, u) + 5ii5jkA{u, t, s) , (A. 8) 

where s, t and u are the usual Mandelstam variables. Note that the phase of eq. ( [A. 8 



corresponds to simply taking the matrix element of the Lagrangian density of a four point 
contact interaction. Projecting out amplitudes of definite isospin yields: 

(s, t, u) = 3A{s, t, u) + A{t, s, u) + A{u, t, s) , 
Tl^{s,t,u) = A{t, s,u) — A{u,t, s) , 

Tl^{s,t,u) = A{t,s,u) + A{u,t,s) . (A.9) 
The needed I = rni ^ KK amplitude can be gotten as: 

T°2(s,t,M) = -VQA{Ti\p^)Tx\p2).K^{vz)K-{vA)) ■ (A.IO) 
We then define the partial wave isospin amplitudes according to the following formula: 

TLA^) = IVf^b £^dcosdPi{cosd)T^,{s,t,u) , (A.ll) 

where 6 is the scattering angle and 



5* is a symmetry factor which is 2 for identical particles (tttt case) and 1 for distinguishable 
particles {KK case). 
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Appendix B 
Chiral Lagrangian 

In the low energy physics of hadrons, it is important to take account of the spontaneous 
chiral symmetry breaking structure. We start here with the U(3)lxU(3)r / U(3)v non-linear 
realization of chiral symmetry. The basic quantity is a 3 x 3 matrix U, which transforms as 

U^Ui^Uul, (B.l) 

where [/l,r ^ U(3)l j^. This U is parameterized by the pseudoscalar as 

U = ^\ e = e'''^/^^ , (B.2) 

where is a pion decay constant. Under the chiral transformation eq. (p.l| ), ^ transforms 
non- linearly: 

e ^ ?7l e i^t (0, Ui^^ Un) = K{^, Uu Un) e ul . (B.3) 



The vector meson nonet is introduced as a gauge field |T7[ which transforms as 



Kp^K^ + '-Kd„K^ , (B.4) 



where ^ is a gauge coupling constant. (For an alternative approach see, for a review, Ref. |T 
It is convenient to define 

which transform as 

Kv^K^ + iKd^K^ . (B.6) 

Using the above quantities we construct the chiral Lagrangian including both pseudoscalar 
and vector mesons: 



^ = -^"^^Tr 



{gp, - v,)'\ - ^Tr [p,p,] - -Tr [F^FM] , (B.7) 
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. , 1. 



where F^^ = dfj,pu — dup^ — ig[Pf_i, Pu] is a gauge field strength of vector mesons. 

In the real world chiral symmetry is explicitly broken by the quark mass term —fnqAiq, 
where m = {rriu + rrid) /2, and Ai is the dimension-less matrix: 



M 



V 



X 



(B.8) 



Here x and y are the quark mass ratios: 



X 



fn 



1 / nid - rriu 

2 V m 



(B.9) 



These quark masses lead to mass terms for pseudoscalar mesons. Moreover, in considering 
the processes related to the kaon, (in this paper we will consider tttt KK scattering 
amplitude.) we need to take account of the large splitting of the s quark mass from the u 
and d quark masses. These effects are included as SU(3) symmetry breaking terms in the 



above Lagrangian, which are summarized, for example, in Refs. Here we write the 

lowest order pseudoscalar mass term only: 



0— mass 



5'Tr 



(B.IO) 



where 5' is an arbitrary constant. 

We next introduce higher resonances into our Lagrangian. First, we write the interaction 
between the scalar nonet field S and pseudoscalar mesons. Under the chiral transformation, 
this S transforms as S* ^ A possible form which includes the minimum number of 

derivatives is proportional to Tr [Sp^pp\ . The coupling of a physical isosinglet field to two 
pions is then described by 

C„ = d^T: ■ df,Ti . (B.ll) 

Here we should note that the chiral symmetry requires derivative-type interactions between 
scalar fields and pseudoscalar mesons. Second, we represent the tensor nonet field by T^y 
(satisfying T^y = Ty^, and T^^ = 0.), which transforms as T^y KT^yK^ . The interaction 
term is given by 

Ct = -72F'Tr [T^yp^py] . (B.12) 

The heavier vector resonances such as p(1450) can be introduced in the same way as p in 
eq. Q. 
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Appendix C 



Unregularized amplitudes 
Amplitudes for the tttt tttt channel 



The current algebra contribution to A{s,t,u) is 



F2 



(C.l) 



The amphtude for the vectors can be expressed in the following form 



Ap{s,t,u) 



2m'i 



t{u ~ -5) _|_ u{t — s 



rni — t m?p — u 



where gpT^T, is the coupling of the vector to two pions. 
For the scalar particle we deduce 

7o - 2^2)' 



Af^{s,t,u) 



2 m%-s 



To calculate the tensor exchange diagram we need the spin 2 propagator 



(C.2) 



(C.3) 



^piUi(^P2V2 ~^ ^piU2^P2'^l) 2^All/^2^l'l!^2 



where 



A straightforward computation then yields the /2 contribution to the tttt scattering 
tude: 



(C.4) 

(C.5) 
ampli- 



Af^{s,t,u) 



16 4 10 2 -'- 2 , / .2 I 2\ 



2fm2 -s) V 3 "^^ ' 3 "^"^ 3' ' 2 



2 vnLs^ 



+ 



3 6mj2 Sm^^^ 



(C.6) 



Amplitudes for tt tt K^K 



Current algebra amplitude: 



2F2 



(C.7) 
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Vector meson contribution: 



Ayector-(7r°7r°,X+X-) 



9k*Ktt 



t{s ~ u) ^ u{s — t) 



mj^, — t mj^» — u 



2\2 



1 1 

+ 



m\* — t 171%;* — u ^ 



Direct channel contribution for the scalar: 



{s — 2m^)(s — 2m 



Cross channel contribution for the scalar: 



Ik* Kit 



{m\ + ml- tf ^ (jnjc + ml- uf 



rriir* — t 

^0 



m\* — u 



Direct channel tensor contribution: 



l2-KTil2KK 



2(m 



h 



Am% 2 2 ; 



+ 



u ml + /;/!■, 



K^m% 2 
2 / 



3 \4m% 2 



2 o +"^^ 



Cross channel tensor contribution: 



X 



+ 



{2ml - s) - 



{2ml - 



2m 



{ml -m\ + 



2ml-. 



-{ml-ml + tf 



{u-ml- ml) + - i'^K - mlf) 



2 
3 



{t-ml- ml) 



2m 



/j.2 / 2 2 

-{t -{m^-m^ 



2\2\ 



+ {t < > u) . 
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